Submicron plasticity: yield stress, dislocation avalanches, and velocity distribution 
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The existence of a well defined yield stress, where a macroscopic piece of crystal begins to plas- 
tically flow, has been one of the basic observations of materials science. In contrast to macroscopic 
samples, in micro- and nanocrystals the strain accumulates in distinct, unpredictable bursts, which 
makes controlled plastic forming rather difhcult. Here we study by simulation, in two and three 
dimensions, plastic deformation of submicron objects under increasing stress. We show that, while 
the stress-strain relation of individual samples exhibits jumps, its average and mean deviation still 
specify a well-defined critical stress, which we identify with the jamming- flowing transition. The 
statistical background of this phenomenon is analyzed through the velocity distribution of short 
dislocation segments, revealing a universal cubic decay and an appearance of a shoulder due to dis- 
location avalanches. Our results can help to understand the jamming-flowing transition exhibited 
by a series of various physical systems. 

PACS numbers: 64.70.Pf, 61.20.Lc, 81.05.Kf, 61.72.Bb 



Understanding the nature of irreversible plastic defor- 
mation is a crucially important issue in current research 
in materials sciences. It has been known for a long time 
that macroscopic size materials begin to yield at a certain 
stress level, depending on several material parameters. 
On macroscopic scales the flowing regime is traditionally 
described by constitutive laws, envisaging plasticity as 
a smooth and steady flow both in time and space. In 
the last decade, however, a completely new picture has 
emerged. By analyzing the emitted sound waves dur- 
ing deformation it was observed that plastic deformation 
is characterized by intermittent bursts of activity [ll-Q- 
Furthermore, recent compression tests carried out on Ni 
microcrystals 0, @] revealed that, if the specimen size 
is in the range of several //m and below, discontinuous 
deformation results in steps of irregular spacing on the 
stress-strain curve. The related strain bursts correspond 
to the sudden, collective motion of line-like lattice de- 
fects, the dislocations. The size distribution of the dis- 
location avalanches decreases by a universal power-law 
with exponent 1.5. The distribution was recovered by 
experiments [1, 0| , by two- (2D) and three-dimensional 
3D) discrete dislocation dynamics (DDD) simulations 
I , and by analytical modeling [l^ . These findings 
indicate that intermittent dislocation avalanches are an 
intrinsic feature of the plasticity of crystals, and their 
size distribution is not affected by the details of the de- 
formation. 

In nanoscale applications this behavior has two im- 
portant consequences. Firstly, the increased relative size 
of the fluctuations makes difficult to control the plastic 
forming process [1] . Secondly, at small specimen sizes the 
yield stress is not well-defined any more. According to 
the conventional definition, the yield stress of a specimen 
is the external stress at 0.2% plastic strain. If the speci- 
men size is below several /im, then, as a result of strain 



fluctuations and statistical effects in the microstructure, 
this value varies specimen by specimen prohibiting 
a material-specific yield stress definition. One can thus 
raise the question, whether it is possible to give a new 
physical definition for the yield stress, or in other words, 
whether the strength of a micron size specimen can be 
defined at all. 

Similar phenomena are observed in completely differ- 
ent physical systems, too. Granular materials, such as 
sand piles, also exhibit yield stress, and at higher external 
driving forces deformation occurs in distinct avalanches 
11 1 . Plate tectonics [T^], fracture dynamics [isj . vortex 
lattices in superconducting films [1J| , and dynamics of 
domains in ferromagnets jl5| are further examples of such 
processes. The common characteristics of these systems 
are marginal stability, power-law distributions without 
characteristic length- or timescales, and driving forces 
that vary much slower than the internal relaxation pro- 
cesses [l6|. For such problems the term "self-organized 
criticality" is widely used 171. Il8|. 



Although irregular plastic response of submicron crys- 
talline materials is by now conceived as a self-organized 
critical phenomenon ]^, [^, [igI . [19| , a physical understand- 
ing, that is, a phenomenology for plastic flow based on 
statistical properties of strain avalanches is still lacking. 
In this paper we present a statistical analysis of the fluc- 
tuating stress-strain response of individual specimens, 
and of the velocity distribution P{v) of short dislocation 
segments. Our main proposition is that the P{v) holds 
the key to many empirical aspects of the flow. It exhibits 
a specific transition at the onset of material yielding seen 
on the average stress-strain characteristics. 

Concerning the dynamics of dislocations, it is com- 
monly assumed that the motion of dislocations is es- 
sentially overdamped. Thus the glide velocity is pro- 
portional to the glide component of the acting Peach- 
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Koehler force [20] per unit length Fg, that is, v = B~^F^, 
where B is the drag coefficient of the dislocations [20j. 
Here Fg — feris + F^, where b is the magnitude of the 
Burgers vector and t^s is the sum of the external stress 
and local resolved shear stress, generated by the other 
dislocation segments via long-range interactions, while 
Fs denotes the short-range forces related to dislocation 
self- interaction and junction formation (20l |. 

To study the submicron plastic response by DDD we 
simulated in 3D several statistically equivalent realiza- 
tions of the yielding of an L = 0.46 ^m edge size alu- 
minium cube subject to compression with constant stress 
rate and free side boundaries (for simulation details see 



2l|). The initial dislocation density was 8 x 10 m 



Part of a typical configuration is shown in Fig. [IJa) [26j. 
Motivated by the fact that in 2D the avalanche size distri- 
bution was found to be very similar to the 3D case 0, 
simulations on 2D systems consisting of straight parallel 
edge dislocations, oriented for single slip with periodic 
boundary conditions, were also performed. A typical 2D 
system is seen in Fig.[TIb) [26]. Although several 3D in- 
gredients, like dislocation multiplication, junctions and 
forest dislocations are absent from the 2D system, there 
is ground to expect that the essential feature governing 
strain avalanches are included in the 2D case. Therefore, 
beside the much less computational cost, an important 
advantage of the 2D simulations is that the role of long- 
range interactions in irregular plastic yielding is isolated, 
so its effects can be better studied. 



(b) 





Figure 1: Dislocation configurations obtained by 2D and 3D 
simulations, (a) Snapshot of a three-dimensional simulation. 
Different colors indicate dislocations on the 12 different slip 
systems, (b) Snapshot of a two-dimensional system consisting 
of straight parallel edge dislocations. 



The plastic strain responses measured during individ- 
ual stress-controlled 3D simulations are seen in Fig. HJa) 
(thin lines). The obtained stress (roxt) vs. plastic strain 
(7) curves exhibit random steps, just like the ones mea- 
sured experimentally by microcrystal deformation 0, Q . 
The plateaus clearly indicate strain avalanches, result- 
ing in different patterns for different samples excluding 
any practical definition of a threshold value. If, how- 
ever, we average over the more than 100 independent 
realizations of submicron samples, the cavalcade of ran- 
dom staircases smoothens into a continuous curve, the 



thick line of Fig. Ufa) Moreover, as seen in the inset of 
Fig. Ufa), there is a threshold stress value Tc « 65 ± 10 
MPa marking the end of the pure power region. The on- 
set of the flow is even more evident in the deformation 
rate 7 vs. external stress relation, plotted in Fig. Hfb), 
undergoing a quite sharp transition at the same Tc- Fur- 
thermore, the root mean square fluctuations over sam- 
ples of plastic strain values at given stresses, depicted in 
Fig.EJc), sharply increase for Text > Tc- 
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Figure 2: Average properties of strain responses obtained by 
DDD simulations in tree dimensions, (a) Stress (rext) ver- 
sus plastic strain (7) curves of individual simulations (thin 
lines). The thick line is the stress-strain curve averaged over 
samples. The black straight line marks the suggested criti- 
cal yield stress. Inset shows the averaged curve on a log-log 
plot, (b) Deformation rate versus external stress, (c) Root 
mean square of the plastic strain fluctuations as function of 
the external stress. 

We now turn to the detailed presentation of the 2D 
results. The simulations were started out from differ- 
ent random distributions with zero net Burgers vector. 
For each of the different 5000 realizations, after letting 
the system relax at zero external stress, stress-controlled 
loading was applied. In contrast to 3D, in the 2D sim- 
ulations all the material parameters can be scaled out 
by introducing natural units as vq = ^/pB~^Gh^ for ve- 
locity, 7o = fo-^p for plastic strain and tq — for 
stress, where p is the total dislocation density and G is a 
combination of elastic moduli [22[. Just like in three di- 
mensions, the individual curves are step- like (Fig. [21[a)), 
but again, a relatively sharp critical transition point can 
be identified on the external stress dependence of the 
average plastic strain (inset of Fig. Ela)), on the defor- 
mation rate 7 (Fig. ^ih)) and also on the plastic strain 
fluctuation (Fig. [3Kc)). For 128 dislocations the critical 
stress obtained is Tc ~ (0.17 ± 0.02) tq. The outstanding 
similarity between the strain responses of 2D and 3D sys- 
tems observed in Figs. [5] and [3] suggests that the simpli- 
fled 2D case, only containing the long-range interactions 
between dislocations, is able to capture the general fea- 
tures of submicron plastic flow. We can conclude that not 
only there is in the statistical sense a smooth stress-strain 
curve for submicron sizes, but also that several corrob- 
orating indicators show the existence of a quite sharply 
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defined threshold stress, presumably a material charac- 
teristics for a given specimen size and initial dislocation 
density. 
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Figure 3: Average properties of strain responses obtained by 
DDD simulations in two dimensions. (a,b,c) Each panel is the 
counterpart of Figs. [5]^a-c) , now obtained by 2D simulations. 
Here the quantities are in natural units defined in the text. 

More can be learned from the detailed analysis of the 
dislocation segment velocity distribution P(v). Since the 
in situ direct measurement of the velocities of large num- 
ber of dislocations is rather difficult, there are no exper- 
imental data available on them so far. Here we report 
about results obtained by DDD simulations in 2D. In 
each individual simulation we determine the distribution 
of the absolute velocities of the dislocation segments, at 
different load levels and take their averaged over 5000 
statistically equivalent realizations (see Fig. HJa)). The 
first thing that needs to be mentioned is the remarkable 
fact that the tails of the distributions exhibit a scale-free 

P{v) « Av-^ (1) 

power-law asymptotic decay with exponent A ~ 3 ±0.02. 




Figure 4: (a) Dislocation velocity distributions at different 
external stress levels. The arrow indicates the direction of 
increasing external stresses, (b) The prefactor yl(rext) (see 
equation (1)) as a function of the external stress, (c) Above 
the yielding transition the velocity spectrum (red) can be de- 
composed into the contribution of avalanches (green) and that 
of the quiescent configurations (blue). 



The theoretical explanation of the cubic decay can be 
obtained under some natural assumptions about the cor- 
relation of dislocations. Previously it was shown that 
the distribution of internal elastic stresses at a random 
location has a cubic tail [1^, presently however we are 
interested in the stress (proportional to the velocity) at 
the positions of dislocations. We give a simplified deriva- 
tion first and then discuss its adaptation to our case. The 
cumulant generating function ^'(g) of the stress distribu- 
tion is given by 

e*(9) ^ J dre*«^p(r) = (Je*«^"=i , (2) 

where the average is taken over the normalized joint dis- 
tribution of iV -I- 1 dislocations /(t-q, . . . tat) and T(r) = 
cos (fi cos 2(p is the stress field generated by a disloca- 
tion, in appropriate units, using polar coordinates. The 
large-r behavior is determined by small g's, so we keep 
the leading term in the Mayer cluster expansion 

^{q) ^nJ (e^'-W - l) /(r|0) dV, (3) 

where the conditional distribution /(r|0) appears. For 
r < qe, where e is small but fixed, the phase factor oscil- 
lates fast so the integral gives at most an order q^ con- 
tribution, if / is nonsingular. Setting the average of r to 
zero, the leading term comes from expanding in q as 

NT A' 
*(9)«-y9'y r^{r)f{r\0)d\^~-q^lnq, (4) 

r>qe 

where A' ~ N J^\os^ tp cos^ 2ip f{r = 0, ip\0) dip , pro- 
vided the small-distance limit of /(r, Lp\0) depends only 
on the angle. Note that 1/f is of the order of area, so 
Nf is nonextensive. As the final step in the deriva- 
tion one straightforwardly shows that the asymptote 
P(r) « A't~^ yields by ^ the g-dependence in (g]). We 
mention that the cubic decay implies a variance logarith- 
mically diverging in the upper cutoff, corresponding to 
the logarithmic factor inffl, which is a typical feature of 
X-ray line-profile tails |24| . 

Now we need to refine the above picture by the fol- 
lowing considerations. Firstly, there are two types of 
dislocations with ± Burgers vectors, so pair correlations 
between all types have to be introduced, as it was worked 
out in [2^. Secondly, we should realize that the equilib- 
rium correlation /eq(r|0) may diverge for small distances 
[2^. However, since dislocations do not move, the stress 
distribution is the Dirac delta centered at the origin, so 
^'oq(<7) = 0. Thus we can use in ^ the time dependent 
deviation from the equilibrium correlation. A detailed 
analysis will be published elsewhere. In conclusion, if 
the time-decaying term in the correlation produces a fi- 
nite amplitude in ^ then the stress and thus the velocity 
distribution must have a reciprocal cubic decay. 
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Returning to the discussion of the simulations, the 
prefactor A in Eq. ([IJ was found to depend on the mo- 
mentary external stress Text- According to Fig. Hl[b), 
^(Tcxt) also undergoes a transition. Although the turn is 
not very sharp, the asymptotes apparently have different 
slopes, whose intersection point again gives Tc of Fig.[3][a- 
c). In a statistical sense, therefore. Text > Tc corresponds 
to the flowing regime, so, for a given size on the average, 
a dynamical transition takes place at Tc, which can be 
considered as a measure of the strength of the material. 
It should be stressed, however, that in an individual sam- 
ple considerable dislocation motion may occur below Tc- 
One can expect that, for samples with macroscopic sizes, 
Tc goes over to the conventional yield stress. Finally we 
note that according to our preliminary investigations the 
rate dependence seems to be weak. 

The investigations on the velocity distributions have 
been repeated also in 3D. Although in this case only a 
smaller ensemble can be afforded and the inherent nu- 
merical noise is much higher, there is a striking similar- 
ity between the 3D and the 2D results discussed above. 
The tail of the distribution is cubic, and its prefactor in- 
creases with the growing external stress. The emergence 
of the cubic decay indicates that its derivation in the 2D 
single-slip case actually has a much broader validity. Fur- 
thermore, as in 2D, the different distributions separate 
after a well-defined point and the small-t! part is nearly 
unaffected by the external stress. We must conclude that 
the observed velocity distribution and its evolution are of 
universal nature. 

The next issue to be considered is whether there is any 
mark of the avalanche activity on the velocity distribu- 
tion of dislocations. In an individual simulation run the 
system alternates between the quiescent and avalanche 
states. This is manifested in surges in the average abso- 
lute dislocation velocity within a sample as function of 
time, the heights identified as avalanches separated by 
low activity periods (see Fig. [5]). These two states can 
be distinguished by thresholding the average dislocation 
velocity in each separate run. At a given external stress 
level it is instructive to calculate separately the averaged 
velocity distributions corresponding to systems being in 
quiescent and avalanche states. The two distributions are 
plotted in Fig. \^c) with blue and green curves, respec- 
tively. As seen, the tail of the total velocity distribution 
comes from systems being in avalanche state. Samples 
in quiescent states contribute only to the small-w part of 
P{v). The superposition of the two convex curves results 
in the characteristic "shoulder" in the distribution func- 
tions plotted in Fig. Hlja). These observations apply in 
three dimensions, too. In sum, the velocity distribution 
exhibits universality in more than one respect. First, the 
low-velocity part of the histogram is nearly independent 
of the external load. Second, the exponent of the decay 
at large velocities appears to be constantly 3 all along the 
loading scenario. Third, the critical yield stress manifests 




Figure 5: Identification of avalanches in individual simulation 
runs. The average of the absolute velocity fluctuates during 
a simulation. In the figure the two curves correspond to two 
different realizations, but all simulation parameters are the 
same. The dashed line denotes the threshold value used for 
avalanche identification. The regions, where an individual run 
is in avalanche state are marked by the corresponding light 
colors. 



itself by an upturn in the amplitude of the power tail, due 
to increased avalanche weight. 

While so far we analyzed peculiarities of the veloc- 
ity distribution of dislocations, Orowan's well known law 
[20j states that the plastic strain rate is proportional to 
the average dislocation velocity, weighted by the Burgers 
vector. So, it is natural to ask what the contribution is 
of the distribution tail obtained above, corresponding to 
avalanches, to the average plastic rate. We found that 
about 80% of the plastic strain comes from the region at 
and above the "shoulder" . At the same time 20% comes 
from the power-decaying part, significant for a tail. This 
implies that the fast dislocations forming avalanches play 
dominant role in the plastic response. 

In summary we should emphasize that small scale dis- 
location simulations abound from the past thirty years. 
On the other hand, experiments were available only on 
macroscopic plasticity, in sizes never reached by simula- 
tions. In the present paper the two approaches meet, now 
that experiments reached down to submicron level, simu- 
lations become more faithful. Increased computer power 
made us possible to describe ensembles never considered 
before, thus a statistical analysis with new conclusions 
and definition of the yield stress in small scale specimens 
could be reached. 
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